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Abstract 



The scaling of pion multiplicities with system size in proton-nucleus and nucleus- 
nucleus reactions at 200 A GeV is investigated. A nuclear geometry calculation 
of incoherent multiple scattering including effects of deceleration of the partici- 
pating particles is presented. Pion multiplicities at high transverse momenta agree 
relatively well with an overall scaling with the number of binary collisions. The 
deviations from this gross scaling for p+A and S+Au reactions at 200AGeV are 
\q ■ compared to model predictions and the implications are discussed. Within this pic- 

■ ture no universal scaling law for high px particle production is obtained. In addition 

to the conceptual problems of such an incoherent model, it is shown that the im- 
plied kinematic relation between a possible p^-broadening and deceleration leads to 
contradictions. The scaling of pion production at intermediate px with system size 
does not indicate a multiple scattering enhancement. 



1 Introduction 



Hadron production at high transverse momenta in very high energy collisions 
is thought to be described by perturbative QCD. The transverse momentum 
and energy dependence of the cross section is however not consistent with 
simple qq scattering processes. Higher order contributions and scaling viola- 
tions have to be included. At somewhat lower energies (y/s < 20 GeV) it is still 
unclear whether meson degrees of freedom contribute to high px scattering via 
the constituent interchange mechanism (for reviews see e.g. [1,2]). 

The interest in high px production has been renewed by the observation of the 
so-called Cronin-effect in p+A collisions at Fermilab [3]. Here it was seen that 
the target mass dependence of particle production when parameterized as A a 
reveals a power larger than one for px > 2 GeV/c, an effect, which was later 
interpreted as the result of multiple scattering of the incident partons [4,5]. 
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In nucleus-nucleus collisions there is still a strong enhancement of particle 
production at high p?- It has been proposed to parameterize the nuclear de- 
pendence as [6]: 



For minimum bias nucleus-nucleus collisions this parameterization yields ex- 
ponents at high pt which are comparable to those in p+A, however, for heavy 
nuclei the statistical significance was so far limited. Recently high precision 
neutral pion data from 200AGeV S+Au collisions have become available [7], 
which corroborate the above result. It is also emphasized in [7] that this en- 
hancement must not necessarily be an isolated feature at high p T , but that 
differences of the shape of the momentum spectra over the whole px range are 
important. 

Still the similarity of the exponents may be accidental, as a lot of information 
from the heavy ion reaction is washed out due to the impact parameter aver- 
aging. In a naive picture, the total cross section in nucleus-nucleus collisions 
should scale like (Ra + Rb) 2 oc (A 1 ^ + B 1 / 3 ) 2 , while the yield of produced 
particles for a given nuclear collision should depend mainly on the apparent 
thicknesses like A 1 ^ ■ B 1 ^'. For proton-nucleus this reduces to a simple depen- 
dence, because A <C B. This is not generally true for a heavy-ion collisions 
(A ^ B), with an exception of the special case where every possible binary col- 
lision contributes the same amount to the total cross section without further 
modification. Then the cross section must be ABa pp .ln case of an anomalous 
nuclear enhancement as observed in [3] such a simple relation does not hold. 

I will try to establish a way of investigating scaling properties which can be 
applied to centrality selected heavy-ion data. Special attention will be given to 
high pt production. There it is reasonable to assume a scaling with the number 
of binary collisions as a baseline estimate. Modifications due to initial state 
multiple scattering should scale with the thickness of the colliding nuclei. The 
rms-thickness properly averaged over all possible nucleon-nucleon collisions 
will be used here. If the nuclear enhancement is caused by initial state multiple 
scattering, one expects the particle yield per binary collision to be a continuous 
and monotonous function of this thickness parameter. 

For the intermediate pt range there have been attempts [8] to describe the 
change in slope of the spectra for heavier systems as a p^-broadening due to 
initial state multiple scattering. There only the shapes of the spectra have 
been compared, no attempt was made to predict the absolute scale of particle 
production. As such a model must also account for changes in the multiplicity 
of particles, it is of interest to simultaneously compare the concepts of multiple 
scattering also to data at lower p T . 
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In this paper I use a simple simulation of incoherent multiple scattering in the 
initial state to predict possible scaling laws for particle production. One has 
to note, however, that the assumption of incoherence and the decomposition 
of the reaction into subsequent binary collisions is very naive and most likely 
not fulfilled. 

I will calculate a certain transverse rapidity kick from the effective thickness 
with the help of a hydrodynamically inspired formula. This recipe may be 
questionable, but is used merely to study the change in shape of the spectra 
- which it should allow - and does not imply any physical significance of the 
hydrodynamic origin of the formula. In addition, a possible energy degradation 
is calculated, and the altered \fs used in the following collisions. This relies 
heavily on the incoherence of the collisions and is in this respect very similar 
to a naive cascade model. It is actually not expected to be applicable in high 
energy reactions. 

Also, the transverse rapidity kick is only relevant for the produced particles. 
The number of collisions the nucleons suffer is calculated by assuming that 
they continue to travel on a straight line. It is further assumed that both the 
rapidity kick and the longitudinal rapidity loss are independent of the y/s of 
the collision in question. This will obviously not be true once y/s gets very 
small. This will be neglected because the concept of energy degradation has 
its problems, and we will only make limited use of it in this paper. For similar 
reasons the effects of varying "source rapidity" depending on the amount of 
multiple scattering will be neglected. 

The predictions will be compared to pion multiplicities in proton-nucleus col- 
lisions at 200 GeV measured at a pseudorapidity rj = 3.26 [3] and S+Au 
data for different impact parameters at 200 A GeV in the rapidity range 
2.1 < y < 2.9 [7]. Such comparisons might help in understanding the limita- 
tions of incoherent ( "random- walk" ) multiple scattering for the description of 
particle production in nucleus-nucleus collisions. 



2 The Model 



The model calculations use the following recipe: For a given reaction system 
nucleons are randomly distributed in the nuclei according to Woods-Saxon 
distributions: 

p(r) oc 



l + exp[(r-r )/0.54] 
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using r = 1.16A 1 / 3 — 1.35v4 -1 / 3 . It is taken care that the hard cores of the 
nucleons (Rh c = 0.6 fm) do not overlap. The two nuclei collide at a given 
impact parameter. A projectile nucleon is supposed to interact with a target 
nucleon, if their minimum distance on straight line trajectories is less than 
their interaction distance (di n t = 1.025 fm).0 The calculations are repeated 
for a large number of random configurations of nucleons within the nuclei. 

From this model the following variables are calculated: 

• The average number of binary collisions N co u, which may happen in a nu- 
clear reaction. This is the main variable controlling particle multiplicities at 
high momenta in a model of independent production. 

• The reduction in particle multiplicity due to the deceleration of participant 
nucleons. It is assumed that a particle suffers a constant rapidity shift Sy 
per collision. The number of collisions a nucleon may undergo before a given 
collision which is responsible for the production of the particle of interest is 
calculated. I will call this the number of prescatterings N pre . It determines 
the total rapidity loss and thereby the center-of-mass energy of the collision 
in question which will in general be lower than the of the incident 
nucleons. From this I can deduce an attenuation factor X in different ways: 

(1) For high pt production it is assumed that the multiplicities scale according 
to (1 — x T ) 9 [9] with x T = 2p T /y/s, so one gets: 
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S-high - ~ Tg- \*) 

(l - 2p T /^) 

(2) For low pt production a scaling with ln(s) is assumed [10]: 

_ 0.88 + 0.44 ln(si) + 0.118 ln 2 ( gl ) 
low ~ 0.88 + 0.44 ln(s ) + 0.118 ln 2 (s ) ' 

(3) The most radical assumption is that of complete absorption. A nucleon 
produces particles only in the very first collision: 

X ahs = { ' Nfn = °' (4) 

Npre ^ 0. 

This corresponds to the limit of maximum momentum transfer in cases 1) 
and 2). 



2 Mostly nucleon-nucleon interactions are presumed, but it should not really matter 
whether the nucleon or the parton picture determines the underlying dynamics. One 
of the assumptions of the model is, however, that all individual binary collisions in 
the nuclei are incoherent. 
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These factors are averaged over all possible binary collisions and attenuation 
factors Xi = (Xi) are obtained. 
• The effective number of prescatterings N e ff. For this, the squareroot of the 
number of prescatterings in each individual binary collision is calculated. 
This is then averaged over all possible binary collisions. It should control 
the effects of initial state multiple scattering. The value of N e ff for a specific 
reaction depends on the way the averaging over the different binary colli- 
sions is done. For each binary collision the corresponding attenuation factor 
Xi has to be taken into account as a weight, because it reduces the multi- 
plicity. This introduces an implicit dependence on the longitudinal rapidity 
shift 5y - a larger deceleration will result in a smaller effective number of 
prescatterings. 

A certain transverse rapidity shift bp is attributed to every prescattering 
- this is the shift apparent in the final particle production, so the rapidity 
shift a single nucleon gets is actually 25 p. Since N e ff already includes the 
effects of random relative orientation of these prescatterings by adding the 
contributions quadratically, the average transverse rapidity shift in a given 
nuclear reaction is just p(N e ff) = 5p ■ iV e //. The effect of the additional 
transverse velocity is then calculated with the help of a formula frequently 
used to describe hydrodynamical transverse flow [11,12]: 

1/pr ■ dN/dp T _ ( p T sinh [p(N eff )} \ ( m T cosh [p(N eff )} \ 

' mT ' ° ^ T ) l \ T ) 

With Ni being the number of pre-scatterings for a given binary collision, 
the two parameters explained above are different averages N pre = (Ni)u and 

N e ff = (^^/N^ \f\ They both have a well defined role within the model studied 

here, however, they also provide in general (esp. in the case Sy = 0) measures 
of the effective thickness of the two colliding nuclei. It is therefore illustrative 
to study particle production as a function of these variables. 

Figure 1 shows sample results of this calculation for p+A and S+Au collisions 
at 20074GeV. The attenuation factor (left side of figure 1) must obviously be 
x = 1 for p+p collisions (i.e. N co u = 1). It decreases for larger system size. 
There is however no unique dependence on N co u, the behaviour is slightly dif- 
ferent for S+Au compared to p+A. One also notices immediately that the 
effects are much stronger for high than for low p? - this is a direct conse- 
quence of the stronger sfs dependence of the high pt production. 

The effective number of prescatterings (right side of figure 1) is also by def- 
inition N e ff = for p+p collisions and increases with system size. There is 
a similar discontinuity when going from p+A to S+Au collisions. The largest 



One should note, however, that while N e jf contains an influence from energy 
degradation there is no such effect included in N pre . 



values for N e ff are obtained without any deceleration effect (5y = 0), a finite 
rapidity loss per collision of by = 0.5 reduces these numbers slightly for low 
Pt and drastically for high p T . 



3 Analysis and Discussion 

The pion cross sections for pt = 0.77GeV/c and pr = 3.08GeV/c from [3] 
have been normalized with the nuclear inelastic cross section for the corre- 
sponding nuclei [13] to obtain multiplicities/event. Neutral pion equivalent 
multiplicities have been calculated by averaging the results for negative and 
positive pions. 

The neutral pion cross sections for S+Au can be very well described by a fit 
with a power- law like formula [7,14]: 

dp 3 \Pt+PoJ 

From this parameterization the pion multiplicities have been calculated at the 
same pt values as above. The systematic uncertainty of the pion yield given 
in [7] will be used as error. 

The reaction cross sections from [7] for the different centrality classes selected 
by the transverse energy Et are used to calculate impact parameter estimates 
assuming a monotonic increase of Et with b. This may lead to relatively large 
uncertainties for the most peripheral reactions which are affected most strongly 
by the experimental trigger bias. While the results presented below include an 
estimate of the systematic error of the model predictions from variations of 
the model assumptions, this may not completely account for the trigger bias 
in the most peripheral trigger class. 

Scaled pion multiplicities are then calculated for different input parameters of 
the model: 

Ncou dp 6 Xi(6y) 

Figure 2 displays the scaled multiplicities at p T = 3.08GeV/c for p+A and 
S+Au as a function of the effective number of prescatterings N e ff. The upper 
left part shows the data without any correction for deceleration effects. One 
immediately realizes that the multiplicity per collision rises with N e ff (i.e. 
with the target mass) in p+A collisions (open circles) - this is just a different 
representation of the anomalous nuclear enhancement. 
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The p+A data have been fitted with equation 5 according to the following 
procedure: The slope parameter T was tuned to describe p+p at high p? 
- this yields T = 240 MeV. With T fixed the other parameters were then 
fitted to the data for the heavier targets in two different ways: a) keeping the 
normalization fixed to describe p+p and b) leaving the normalization as a free 
parameter. Both curves are shown in figure 2. The fit (a) does not correctly 
describe all the p+A data, while fit (b) accounts well for the heavier targets 
but misses the p+p point. The same analysis can be done assuming a finite 
amount of deceleration of the participant nucleons, the qualitative picture of 
the dependence being unaltered. The fit results are summarized in table 1. 
It can be seen that the necessary transverse rapidity shift Sp increases with 
larger (longitudinal) rapidity loss Sy of the participants. One should note that 
for by — 1 a transverse rapidity shift of Sp ~ 7 would be necessary to account 
for the rise in the data. This is completely unphysical. 

In this context it is of interest to compare the parameters necessary to ac- 
count for the enhancement with the kinematics of scattering processes. Every 
gain in net transverse momentum p' T is related to a loss in momentum p' L 
available for particle production. For elastic scattering one can easily deduce 
from momentum conservation a relation between the transverse rapidity shift 
Speiast a nucleon obtains in a collision and the longitudinal rapidity loss Sy e i ast 
it suffers: 

[tanh(?/o)] 2 = [tanh(y - Sy e i ast )] 2 + [tanh((5p e i ast )] 2 . (8) 

This relation is displayed in figure 3 as a solid line. For any inelastic process, 
the transverse rapidity can only be smaller than this limit for a given lon- 
gitudinal rapidity loss. This allowed region is indicated in figure 3 as a grey 
area. As a good approximation this relation can also be used as a limit for the 
transverse shift obtained in a multiple scattering picture.0 The symbols show 
the parameter values extracted from the fits of equation 5 to the data. It is 
obvious that these values are not consistent with scattering kinematics. The 
maximum transverse rapidity shift allowed for a given 5y would not suffice to 
explain the enhancement observed in p+A collisions. 

This observation can most easily be interpreted as a failure of the underly- 
ing picture of incoherent scatterings. This is actually expected because of the 
short time scales involved in hard scattering processes. In a coherent model 
there is no possibility of an energy degradation between two "individual" col- 
lisions like it is implemented in the model used here, and thus the kinematical 
considerations do not apply in the same way. The results obtained may be 

4 The estimate would be different for cases, where the collision axis of the binary 
collision is strongly tilted relative to the overall collision axis - this will not be 
considered here. 
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regarded as a warning against too naive use of such simplified models. I will 
therefore not discuss these cases further. 

The S+Au data have also been included in figure 2. For all cases the data 
agree with the extrapolation from p+A for intermediate centralities, however, 
the increase with N e ff is much stronger than predicted by the extrapolation. 
The difference in system size between p+A and S+Au data may be relevant: 
While the parameter N e ff is comparable e.g. for p+W and peripheral S+Au, 
the number of binary collisions in peripheral S+Au is larger by a factor of 
two. 

Ignoring the data points for the peripheral S+Au collisions the heavy ion data 
would be consistent with the heavier target p+A data. In this case, however, 
the p+p data cannot be described. There is no universal quantitative descrip- 
tion of both p+A and S+Au data within this multiple scattering picture, 
although it seems to be able to predict the scaling behaviour for both groups 
of data qualitatively. 



Still, as calculations using a "random-walk" picture have been advocated for 
the description of nucleus-nucleus collisions [8], I will compare the present 
model to particle production at lower transverse momentum (p T = 0.77 GeV/ c) 
Figure 4 shows the scaled multiplicities for this case. Without any deceleration 
effects (upper left) the multiplicity per collision decreases for increasing target 
size in p+A reactions. Like for the high case this reflects the value of the 
"Cronin" -exponent a - here one finds a < 1, which indicates that the nuclei 
are not at all transparent for these softer processes. 

Introducing a correction for the attenuation with intermediate values of Sy = 
0.25 — 0.5 does not really compensate for the decrease. If one assumes complete 
"absorption" , so that only the first binary collision contributes to particle pro- 
duction, the scaled multiplicities show a small increase The particle yields at 
low pt in p+A can obviously be understood without any strong enhancement 
via initial state multiple scattering. They are very close to the expectation 
from production only at the surface of the nuclei. The multiplicities in S+Au 
collisions show an additional increase compared to the expectation from p+A. 



5 A rapidity loss 5y > 2 is necessary to compensate for the decrease. This is for 
practical purposes very similar to the complete absorption case. 
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4 Conclusions 



Pion production at high px scales approximately with the number of binary 
collisions. The enhancement relative to this scaling in p+A and S+Au has been 
studied as a function of an effective thickness parameter. Results indicate that 
it will be difficult to describe the enhancement in both data sets with a unique 
function of the thickness. 

The description of the nuclear enhancement at high px in a picture of inco- 
herent scatterings ( "random- walk" ) in the present calculation would require 
parameters which contradict momentum conservation in individual scatter- 
ings. This confirms that the p^-broadening at high px from initial state multi- 
ple scattering should not be described incoherently. Recent perturbative QCD 
calculations including the effects of intrinsic pt and initial state multiple scat- 
tering [15] were able to reproduce data for neutral pion production from central 
S+S reactions [7] as well as preliminary results for central Pb+Pb collisions 
[16,17]. It would be of interest to see, whether a QCD calculation with co- 
herent multiple scattering effects would be able to describe also the scaling of 
particle multiplicity with centrality in heavy ion collisions. 

Particle production at intermediate px in p+A is consistent with production 
only at the nuclear surface. In this scenario there is little room for an enhance- 
ment by initial state multiple scattering. It must therefore be considered very 
unlikely that e.g. the increase in inverse slope of the momentum spectra with 
system size also at intermediate px can be attributed to this mechanism in 
the way discussed in [8]. 

A simple explanation of increasing slopes with larger system size might be 
an admixture of the two components, where the soft component dominates at 
low pt with increasing importance of the hard component with increasing px- 
The ratio hard/soft would roughly scale as volume/surface and would thus 
increase the high p T part of the spectra more rapidly with increasing system 
size. This could explain the behaviour in p+A collisions at intermediate px- 
The particle production in S+Au collisions might call for additional thermal 
or hydrodynamical production not present in p+A, while for the explanation 
of the high px production more refined calculations of hard scattering appear 
to be necessary. 
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Table 1 

Transverse rapidity shift per collision dp from a fit of equation 5 to data at px = 
3.08GeV/c for p+p and p+A reactions at 200 GeV [3]. 
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Fig. 1. Parameter values obtained in the model calculation for p+A and centrality se- 
lected S+Au collisions. Left: attenuation factors Xj for low- and high-p^ production 
assuming Sy = 0.5, right: effective number of prescatterings N e ff calculated without 
deceleration (Sy = 0) and for low- and high-p^ production assuming Sy = 0.5. The 
filled symbols show results relevant for high pr, the open symbols for low px- 
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Fig. 2. Normalized particle multiplicities at pr = 3.08GeV/c for p+A and S+Au 
collisions as a function of the effective number of prescatterings N e ff. The lines 
show fits of the multiple scattering model (equation 5) to the p+A data including 
(a) and excluding (b) the p+p data point. 
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Fig. 3. Transverse rapidity 5p as a function of the longitudinal rapidity loss Sy for 
nucleon-nucleon scattering processes at yfs = 19.4 GeV. The solid line shows the 
relation for elastic scattering (equation 8), the symbols show the fit parameters of 
the multiple scattering model given in table 1, where the black squares make use of 
the p+p data point and the grey circles are obtained without it. 
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Fig. 4. Normalized particle multiplicities at px = 0.77GeV/c for p+A and S+Au 
collisions as a function of the effective number of prescatterings N e ff assuming no 
rapidity loss (left) and complete absorption (right). The lines show exponential fits 
to guide the eye. In the representation of the data in the right plot the same values 
as for 5y = have been used for N e ff, because the true values would be N e ff = 0. 
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